Breakdown of weak-field magnetotransport at a metallic quantum critical point 
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We show how the collapse of an energy scale in a quantum critical metal can lead to physics beyond 
the weak-field limit usually used to compute transport quantities. For a density-wave transition we 
show that the presence of a finite magnetic field at the critical point leads to discontinuities in the 
transport coefficients as temperature tends to zero. The origin of these discontinuities lies in the 
breakdown of the weak field Jones-Zener expansion which has previously been used to argue that 
magneto-transport coefficients are continuous at simple quantum critical points. The presence of 
potential scattering and magnetic breakdown rounds the discontinuities over a window determined 
by t A < f where A is the order parameter and r is the quasiparticle elastic lifetime. 

PACS numbers: 75.10.Lp, 72.15.Eb, 71.10.Hf, 71.18.+y 
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Quantum critical points in metallic systems have gen- 
erated considerable theoretical and experimental interest 
in recent years Q. They are realized by tuning the fi- 
nite temperature critical point of a phase transition in a 
metal to zero temperature. A variety of tuning param- 
eters have been used including hydrostatic pressure 0, 
chemical composition Q and quantum criticality can oc- 
cur serendipitously p|. In insulating quantum magnets, 
magnetic fields often provide a "handle" by which they 
may be made quantum critical [f| and there has recently 
been much interest in doing the same in the metallic case. 
Notable examples include the metamagnetic quantum 
critical endpoint seen in Sr3Ru207 and also the an- 
tiferromagnetic quantum critical point in YbRh2Si2 0. 
The collapse of the characteristic energy scale for fluctu- 
ations can induce deviations from Landau Fermi liquid 
theory 0, Q seen, for example, in transport quantities. 
These are usually computed by assuming linear response 
to driving fields, yet the vanishing energy scale could in- 
validate that assumption [Io| . 

In this paper we explicitly demonstrate such a break- 
down by considering magneto-transport in the simplest 
class of quantum critical metal tuned by magnetic field. 
The presence of the small but finite magnetic field at 
the transition point leads to modified transport response. 
This has assumed added significance recently because of 
the suggestion that the discrepancies between theory and 
experiment in quantum critical matter Q may originate 
in part from a transition at the critical point between 
localized and de-localized spins and with a concomitant 
change in the Fermi surface volume It is suggested 
that the Hall coefficient would reveal this volume change 
and recently evidence ^3] f° r such a change has been 
presented at the field driven quantum critical point in 
YbRli2Si2. Here, we consider the T — > limit of a 
spin or charge-density wave transition in finite magnetic 
field, so the quantum critical fluctuations play a small 
role in quasi-particle scattering when compared with elas- 
tic scattering from impurities. This limit has been con- 
sidered by Bazaliy et al. 0, Q] recently for non-field 
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FIG. 1: (a) Summary of magnetotransport at a zero tem- 
perature density wave (DW) transition as a function of the 
dimensionless gap, rA, and magnetic field, u t. In the re- 
laxation time approximation there is a discontinuity in trans- 
port coefficients along the line rA = due to the collapse 
of the weak-field regime. Disorder and magnetic breakdown 
will smooth this discontinuity into a crossover. A field-driven 
quantum critical point is characterize by a trajectory through 
this phase diagram: (b) a magnetic field suppressing the DW 
phase, or (b) stabilizing the DW phase. 



driven quantum critical points under the assumption that 
a weak-field expansion in magnetic field can be made. 
They conclude that there are no anomalies in the Hall 
conductivity and that, in the absence of perfect nesting, 
changes in other transport quantities are generally linear 
in the energy gap so change continuously at the critical 
point. 

The essential result of our work is that at a density 
wave (DW) quantum critical point the weak-field (Jones- 
Zener) expansion — upon which previous magnetotrans- 
port work was based — breaks down because of the van- 
ishing DW gap A. This breakdown leads, in the re- 
laxation time approximation, to discontinuities in the 
magnetoconductivities: jumps in the resistivity and Hall 
coefficients at the transition, despite the fact that the 
phase transition itself is continuous. For magnetic fields, 
B > A/(evpT) we show that the magnetoresistance is 
non-analytic in the field (~ \B\) — a result not obtain- 
able from a weak-field expansion. While our results are 
valid at any zero temperature DW transition, they as- 



sume particular significance at a field driven quantum 
critical point where A is tuned to zero at finite magnetic 
field. The interdependence of the gap and magnetic field 
define a trajectory through Fig. QJa) [shown schemati- 
cally in Figs, [Hb) or (c)] which inevitably explores the 
non-analytic region. We also go beyond the relaxation 
time approximation to consider the effects of disorder 
and magnetic breakdown very close to the critical point. 
While these smooth the discontinuity they do not alter 
the magnitude of the changes in conductivity across the 
transition at a finite magnetic field. 

We consider the following highly simplified model of 
transport near a density- wave transition (DW). We take 
the Fermi surface in the paramagnetic state to be circu- 
lar (and assume two dimensionality for ease of computa- 
tion) . We then consider density wave formation to induce 
a periodic potential on the otherwise free electrons and 
assume this potential to be proportional to the mean- 
field order parameter: V± Q (B) = A Q ^(B C - B)/B c . 
Here, the periodic potential (assumed real so both signs 
of Q are present) has Fourier components at the order- 
ing wave-vector Q of the DW and its perturbation on 
the free electrons describes Bragg scattering off the DW. 
Degenerate perturbation theory in the presence of this 
potential will gap the dispersion whenever a resonance 
condition is met: eg = e%± n Q- Since transport in the 
T — » limit that we will be considering is dominated by 
the Fermi surface, we need only consider the most im- 
portant Bragg scattering matrix element. For example, 
near the point e k = £k-Q the dispersion is modified to be 



E k = \ (ejfc + e k _ Q ) ± y 1 (e k - e fe _ Q ) + V* . (1) 

When the chemical potential falls in one of the gaps, as it 
will if 2kp > Q, then the Fermi surface is modified by the 
density wave and it is the change in transport properties 
induced by this modification that we wish to study. 

In Fig. El we illustrate the change in Fermi surface 
topology that we envisage. There are a number of possi- 
ble scenarios 

• 2kp < Q: no gaps appear at the Fermi surface and 
transport coefficients will change smoothly through 
the critical point. We therefore do not consider this 
case here. 

• 2kp > Q, Q incommensurate: gaps appear on the 
Fermi surface and lead to an open section of Fermi 
surface [Fig. EH]- In the absence of twinning, trans- 
port properties would discriminate between cur- 
rents parallel and perpendicular to Q. 

• 2kp > Q, Q commensurate with the lattice (e.g. an 
antifcrromagnet). Transport is isotropic through 
the transition and the Fermi surface remains closed. 
We treat this case by considering two real DWs 



2 




FIG. 2: Fermi surface reconstructions at a density wave tran- 
sition. The solid lines indicate the Fermi surface to be used for 
computing transport, while the dashed Fermi surfaces show 
the dispersion in an extended zone scheme, (a) An incom- 
mensurate DW transition leads to open Fermi surface sheets, 
(b) At a commensurate DW transition (Q% = Qi + g) the 
Fermi surface can break into closed pockets. The small gap 
near the DW transition gives a local radius of curvature of 
the Fermi surface ~ A/vf which leads to a breakdown of the 
weak- field expansion in the limit A — *• + . 

with Q = (±Q x ,±Q y ) and Q = (±Q x ,TQ y ) (see 

Fig.Eb)- 

• 2k F = Q (nested) [H. Since the perfect nesting 
not be met very close to the critical point Q this 
case will ultimately reduce to one of the cases above 
and is therefore not considered here. 

For the initial analysis we use the classical Boltzmann 
equation in the relaxation time approximation r: we are 
envisaging the T = limit of the conductivities and ig- 
nore inelastic processes. Rather than solve the Boltz- 
mann equation order by order in the magnetic field fl5| 
we solve to all orders in the field directly using the 
Chamber's formula [l7|: 

For each area element of the Fermi surface, dS, we in- 
tegrate the velocity, v(t) measured along a semiclassi- 
cal quasiparticle orbit. These orbits are defined by the 
Lorentz equation of motion 

dk 

h— = -evxB, (3) 
dt 

where v — V k^(k) / h. In this paper we will always assume 
the magnetic field is perpendicular to the 2D electron 
fluid. 

Our main numerical results are illustrated in Fig. 
The circles illustrate the conductivities obtained by in- 
tegrating the Chambers formula when an DW gap with 
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FIG. 3: The components of the conductivity tensor (a) a xx , 
(b) o X y and (c) a vv (in units of crn = ne 2 r/m) computed 
when a mean-field like DW gap is pushed to zero at a mag- 
netic field equivalent to uj t = 0.2. The circles are the numer- 
ical solution of the Chamber's formula within the relaxation 
time approximation (RTA). The dashed lines show the ana- 
lytic results. The solid line goes beyond the relaxation time 
approximation by including the role of disorder in washing 
out the DW gap. 



Q parallel to the y axis, is suppressed to zero at a criti- 
cal field corresponding to uo t = 0.2 with A /ep = 0.01. 
Here cj = eB/m, the precession frequency around the 
whole circular Fermi surface in the absence of a DW and 
the transition is occurring in the region where uj t ^ 1 — 
i.e. naively the weak field regime. We clearly see discon- 
tinuities in the conductivities a xx and <j xy (and a very 
small one in a yy ) at the critical field. In the case of 
closed Fermi surfaces (not shown) a clear discontinuity 
occurs in all three conductivities. Having observed that 
the conductivities can be discontinuous at a DW transi- 
tion in a magnetic field we now consider the underlying 
physics. 

Within Boltzmann transport theory the effect of a 
magnetic field on the conductivity is caused by the pre- 
cession of the out-of-equilibrium distribution around the 
Fermi surface according to Eq. |2| between scattering 
events. Quasi-particles are deflected by the local Hall 
angle which is related to the local radius of curvature of 
the Fermi surface (R) and the Lorentz equation of mo- 
tion: 



5k = ROn = evpBr 



?h = 



evpB 



(4) 



longitudinal conductivity Atr xx /a xx (0) ~ cos8h — 1 ~ 
9 2 H giving the usual magneto conductance quadratic in 
applied field. The Hall angle must be small for weak- 
field response. However at a DW transition a gap opens 
in the Fermi surface leading to a local radius of curvature 
of the Fermi surface of order A/vp which vanishes at the 
quantum critical point itself [see Fig. [He)]. The condi- 
tion for being in the weak-field regime is therefore that 
-Bwcak-ficid < A/(evpr) or equivalently uj q t < rA/(kpl) 
where / is the mean free path. This upper bound vanishes 
at the critical point so magnetotransport there is never 
weak-field. Thus the results of Bazaliy et al. are valid 
only provided the limit B/A — * can be taken which 
is certainly not the case at a field-driven DW quantum 
critical point. 

Instead magneto-transport is dominated by a fraction 
of the quasiparticles (proportional to evpBr) that are de- 
flected around the cusp in the Fermi surface between scat- 
tering events. This leads to a magnetoconductance pro- 
portional to \B\ — i.e. non-analytic in the field and thus 
beyond any weak-field expansion. This physics is beau- 
tifully illustrated by the magnetoresistance of a square 
Fermi surface calculated by Pippard |isj . 

We can use Pippard's method to compute the size of 
the discontinuities in all three independent components 
of the conductivity tensor. In the limit that A — > + 
we can solve the Boltzmann transport equation on each 
segment of the Fermi surface and then combine the solu- 
tions by ensuring that the magnitude of the distribution 
function is continuous. The result is that 
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This deflects the current, leading to a reduction of the 



On the paramagnetic side (A — > 0~) the conductiv- 
ities are as above but without the terms proportional 
to P(cx,u/ t). Thus the all three magnetoconductivities, 
Vxx, &xyi vyy, show discontinuities with magnitude of 
order, lo t, (lo t) 2 and (w t) 3 respectively so the most 
dramatic jump is the magnetoresistance footnote. These 
analytic solutions are shown as the dashed lines in Fig.|3 
We now extend our calculation beyond the relaxation 
time approximation. If tA < 1 this calculation cannot 
be valid since the quasiparticles would be unable to no- 
tice the Bragg scattering from the DW above the scatter- 
ing from impurities. One would expect this to washout 
the discontinuities in magnetoconductivities since quasi- 
particles will tend to remain on original Fermi surface. 
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Magnetic (Zener) breakdown when B > 2A 2 /(evp) 0] 
as the same effect. (See Green and Sondhi llO| for E 
field breakdown.) We include this phcnomcnologically 
in our calculation in the following fashion. Rather than 
hybridize the dispersion (as in Eq. ^) we maintain a cir- 
cular Fermi surface and treat the DW potential in the 
collision integral as a resonant scatterer between points 
of the Fermi surface that satisfy the resonance condition. 
This gives the following transport equation for the out of 
equilibrium distribution function g(e,9) 

dg 



ev- Erdefo + ^oT— 
oO 



g- (tA) 2 [5(cos6 + r]) 
+5{co8 6-r})]\g(e,6)-g{e,ir-6)] , (9) 



where /o is the Fermi function and r\ = Q/2k. Since the 
solution is periodic around the Fermi surface it may be 
solved by Fourier transform and, for a circular Fermi sur- 
face gives the following expression for the conductivities 
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K(a, uo q t) 



2Au> TP(a, uj t) cos 2 a 
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(10) 
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and A 



4(tA) 2 



These conductivities are shown as solid 



lines on Figs. [31 Note how these expressions interpolate 
between the paramagnetic solution when tA = (no DW 
scattering), and the Pippard result in the limit tA » 1 
where the angle dependent scattering effectively mimics 
the reconstructed the Fermi surface. Thus we see that 
disorder washes out the discontinuity over a region in 
field determined by rA < 1. 

We have also considered the case of closed Fermi sur- 
faces. In that case a xx = a yy and both show a disconti- 
nuity of fractional order lu t. The other difference from 
the case of open Fermi surfaces is that, as for all open 
Fermi surfaces, a yy remains finite in the high field limit: 
lu t 3> 1. This is the regime where Landau level quan- 
tization of the closed Fermi surfaces would also become 
important and is not considered here. 

In summary, we have shown that at a simple density 
wave (DW) quantum critical point the weak-field regime 
of magnetotransport collapses to zero field with the size 
of the gap. At finite field in a clean metal one would 
expect to see discontinuities in the magnetoresistance of 
order of the magnitude of the Hall angle: ui q t. This ef- 
fect will be significant at a field driven quantum critical 
point where by definition the field is finite at the transi- 
tion point. The case of YbRh2Si2 already shows features 



in the Hall effect which suggest that it falls outside the 
class of DW quantum critical points. A prediction from 
this work is that non-field driven DW quantum critical 
points should should a low field cross-over to a transverse 
magnetoresistance that is linear in the applied field. It 
would be interesting to look for such an effect in the 
Cr\^ x V x system under pressure [2(| where we estimate 
that rA = 1 at (x — x c )/x c ~ 0.1 and u> t ~ (fci?Z) _1 
at B ~ IT. Very recently we have learnt that G&^Kx^O-j 
shows exactly the linear magnetoresistance we predict 
and is argued to be a small gap density wave system [2lJ . 
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